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Abstract. — Two-dimensional layered aperiodic Ising systems are studied in the extreme 
anisotropic limit where they correspond to quantum Ising chains in a transverse field. The 
modulation ol the couplings follows an aperiodic sequence generated through substitution. Ac- 
cording to Luck's criterion, such a perturbation becomes marginal when the wandering exponent 
of the sequence vanishes. Three marginal sequences are considered: the period-doubling, paper- 
folding and three-folding sequences. They correspond to bulk perturbations for which the critical 
temperature is shifted. The surface magnetization is obtained exactly for the three sequences. 
The scaling dimensions of the local magnetization on both surfaces, Xma and Xms, vary con- 
tinuously with the modulation factor. The low-energy excitations of the quantum chains are 
found to scale as with the size L of the system. This is the behaviour expected for a strongly 
anisotropic system, where z is the ratio of the exponents of the correlation lengths in the two 
directions. The anisotropy exponent z is here simply equal to x,ns +Xm3- The anisotropic scaling 
behaviour is verified numerically for other surface and bulk critical properties as well. 



Resume. — On etudie des systemes d'Ising aperiodiques en couches a deux dimensions dans 
la limite anisotrope extreme oil ils correspondent a des chaines quantiques d'Ising en champ 
transverse. La modulation des interactions est engendree par une suite aperiodique obtenue 
par substitution. D'apres le critere de Luck, une telle perturbation devient marginale lorsque 
I'exposant de "divagation" associe a la suite s'annule. Trois suites marginales sont examinees : 
"doublement de periode", "phage de papier" et "phage ternaire". EUes correspondent a des 
perturbations de volume entrainant un changement de temperature critique. Des expressions 
exactes de I'aimantation de surface sont obtenues pour les trois suites. Les dimensions anormales 
de I'aimantation locale sur les deux surfaces, Xms et Xms, varient continiiment avec le facteur 
de modulation. Les excitations de basse energie des chaines quantiques varient en avec la 
taille L du systeme. C'est la le comportement attendu pour un systeme fortement anisotrope 
oil z est le rapport des exposants associes aux longueurs de correlation dans les deux directions. 
L'exposant d'anisotropie z s'exprime simplement comme la somme des exposants magnetiques 
Xms + Xms- Le comportcmeut d'echelle anisotrope est verifie numeriquement pour d'autres 
proprietes critiques de surface et de volume. 



(*) Author for correspondence (e-mail: turban@lps.u-nancy.fr) 
(**) Unite de Recherche Associee au CNRS (URA 155) 
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1. Introduction. 

The discovery of quasicrystals Q has opened a new area of research which has been quite 
active in the last years (see [2-6] for reviews). In the field of critical phenomena, quasiperiodic 
or aperiodic systems are quite interesting since they offer the possibility to interpolate between 
periodic and random systems. 



The Ising model [7-9], the percolation problem |10, |Tl| and the statistics of self-avoiding 
walks [Q were found to display the same critical behaviour on the two-dimensional {2d) 
Penrose lattice as on a periodic one. Universal behaviour was also obtained in 3c? jl3j. On 
the other hand, an aperiodic modulation of the couplings was shown to influence interface 
roughening in 2d [|l^, |l^: a continuously varying roughness exponent was obtained with the 
Fibonacci sequence. 

Layered systems provide another type of aperiodic system which is now accessible experimen- 
tally in three dimensions, via artificially grown multi-layers. Long after the pioneering work of 
McCoy et al [16-19] on the 2c? randomly layered Ising model, some results have been obtained 
with an aperiodic modulation of the interlayer couplings [20-25]. In most cases the specific 
heat displays the Onsager logarithmic singularity but some aperiodic sequences have also been 
found for which the singularity is washed out [|2^, like in the randomly layered system. 

The situation remained unclear until Luck recently proposed a relevance-irrelevance crite- 
rion [|27|| adapted to aperiodic perturbations. As in the Harris criterion for random systems |2^ , 
the strength of the fluctuations of the couplings, on a scale given by the correlation length, 
is of primary importance for the critical behaviour. Thus an aperiodic perturbation can be 
relevant, marginal or irrelevant, depending on the sign of a crossover exponent involving the 
correlation length exponent of the unperturbed system v and the wandering exponent w which 
governs the fluctuations of the aperiodic sequence [^, ^ . The criterion has been extended to 
the case of strongly anisotropic systems with uniaxial aperiodicity [^ , explaining the interface 
roughening results, and later generalized to c?-dimensional aperiodic systems [|3^ . 

Some exact results for the surface magnetization of the 2c? layered Ising model in the extreme 
anisotropic limit [|33|[ have been obtained with irrelevant, marginal and relevant aperiodic per- 
turbations [34-37]. The critical behaviour is as expected on the basis of the Luck criterion: 
one obtains stretched exponentials or flrst-order surface transition when the perturbation is 
relevant and power laws with continuously varying exponents in the marginal case. Conformal 
aspects have been also discussed psf . 

In the present work, we continue our study of marginal sequences, through a detailed analysis 
of the strongly anisotropic scaling behaviour which has been recently found in the layered Ising 
model [^ . In section 2, we recall some generalities about the layered Ising model in the extreme 
anisotropic limit and aperiodic sequences generated through substitution. The surface and 
bulk properties of the layered Ising model with period-doubling (PD), paper- folding (PF) and 
three-folding (TF) modulations are presented and analysed in the light of anisotropic scaling 
in section 3. Some properties of the PF and TF sequences are discussed in two appendices. 



2. Layered Ising model and substitution sequences. 

2.1 Aperiodic Ising quantum chain. — We consider a 2c? layered Ising model with constant 
interaction Ki along the layers and modulated interactions K2{k) between successive layers, k 
and fc-|-l as shown in figure 1. The extreme anisotropic limit corresponds to Ki — s- oo, K2 0, 
while keeping the ratio Xk = K2{k)/Kl fixed, where = — ^ Intanhi^i is the dual coupling. 
In this limit, the row-to-row transfer operator T — exp[—2Kl7{\ involves the Hamiltonian of a 



Fig. 1. — Layered aperiodic Ising model in the extreme anisotropic limit (top) and the corresponding 
quantum chain (bottom) 



quantum Ising chain in a transverse field: 



2 

A; 



where the tr's are Pauh spin operators and is the inhomogeneous couphng defined above. 

Using the Jordan- Wigner transformation and a canonical transformation to the diagonal 
fermion operators rja, the quantum chain Hamiltonian can be put into diagonal form 

n = Y.eJ^U^-l) . (2) 



2 

The fermion excitations Eq follow from the solution of the system 

eafa(fc) = -<Pa(fc) - Afc<^„(fc+1) 

e„<Z'„(fc) = -Afe-i<Z'„(fc - 1) - tf'„(fc) (3) 

which may be rewritten as a single eigenvalue equation, either for # or The components 
(I^a{k) and ^a{k) satisfy appropriate boundary conditions and the eigenvectors are assumed to 
be normalized. 

With an aperiodic modulation, one may write the couplings as: 

Afe = Xrf" , (4) 
where fk, which may be or 1, is determined by the aperiodic sequence. Let 

L 

riL^^fk (5) 

k=l 

be the number of modified couplings on a chain with length L; their asymptotic density is 
Poo = liniL^oo "-l/-^- The critical coupling Ac is such that limi^oo(l/i) X]fc=i Afe = 0, 
which, using equation (^, gives: 

Ac = r-P- . (6) 
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2.2 Substitution sequences and Luck's criterion. — The sequences considered below 
are generated via substitutions, either on digits or on pairs of digits. 

The PD sequence follows from the substitutions <S(0) = 11 and 5(1) = 10, which, 
starting on 1, give successively: 

n = 1 1 

n = 2 1 

n = 3 10 11 

n=A 10 1110 10 

Details about the properties of the period-doubling sequence can be found in p7|. 



Two other sequences will be considered in the following: the PF sequence WJ] generated 
through substitutions on pairs of digits, which is studied in appendix A, and the TF se- 
quence ]45| ] studied in Appendix B. 

Most of the properties of a sequence can be deduced from its substitution matrix |3^ 
with entries My giving the numbers nf*''''' of digits (or pairs) of type i in S{j). In the case of 
the PD sequence, one obtains: 

M^(|:i:::)^(^!)' 

One may verify that the entries in A£" give the same numbers of digits (or pairs) in the 
sequence after n iterations. Let be the right eigenvectors of the substitution matrix and A^ 
the corresponding eigenvalues. The asymptotic density of digits (or pairs) of type i are related 
to the components of the eigenvector corresponding to the leading eigenvalue Ai through 

allowing the calculation of the critical coupling via equation (^ . The length of the sequence 
after n iterations is related to the leading eigenvalue through i ~ A" so that Ai > 1. The 
cumulated deviation A(i) from the average A at the length scale L reached after n iterations 
may be expressed as 

L 

A(L) = ^(A, - A) = A(r - l){nL - Lp^) ~ ^jAaj" ~ <5L" (10) 
fc=i 

where ^ = A(r — 1) is the amplitude of the modulation and uj is the wandering exponent given 
by: 

ln|A2| 1^^. 

Thus the mean shift of the coupling strength at the length scale L is proportional to L'^^^. 

Near the critical point of the pure system, the relevant length is the correlation length 
^ ~ and the aperiodicity induces a shift in the critical temperature 5t ^ ^"^-i n-j ^-i^(i^-i) 
to be compared with the deviation from the critical point t p7| . One obtains a ratio 

= i + K^-i), (12) 



which is divergent when ^ > 0, which corresponds to a relevant perturbation. When (f) < 0, the 
perturbation, which is washed out when the critical point is approached, is irrelevant. When 



Fig. 2. — Non-universal scaling dimension Xms of the surface magnetization as a function of the mod- 
ulation factor r for the period-doubling (PD) , paper-folding (PF) and three- folding (TF) sequences. 



Working in the basis which diagonalizes the Hamiltonian 



one obtains: 



(13) 



where Eq is the ground-state energy. Rewriting the matrix element (i|cr^|0) and the excited 
state |i} in terms of diagonal fermions leads to: 



One may notice that one-fermion excited states alone contribute to the sum. 

In the ordered phase, A > Ac, the lowest excited state \a) =77! |0) becomes degenerate with 
the ground state and gives the only nonvanishing contribution to when r ^ 00 in It 
follows that TOg is given by the first component ??i(l) of the normalized eigenvector associated 
with the lowest excitation ei 




(14) 



a 
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Table I. — Surface energy exponent for the three marginal sequences studied as function of the 
modulation factor r. For each sequence Xes is the extrapolated numerical value whereas the next 
column gives the conjectured analytical result. The figure in brackets gives the uncertainty on the last 
digit. 





period-doubling 


paper-folding 


three-folding 


r 














0.2 


2.396 71(1) 


2.396 711 


4.009(1) 


4.008 961 


1.962(3) 


1.962 149 


0.5 


2.076 34(4) 


2.076 341 


2.669 92(5) 


2.669 925 


1.882 1(2) 


1.881 995 


1 


2.000 0(4) 


2 


2.000 0(3) 


2 


2.000 1(2) 


2 


1.5 


2.026 3(1) 


2.026 274 


1.766 41(8) 


1.766 412 


2.156 873(3) 


2.156 866 


2.0 


2.076 4(3) 


2.076 341 


1.669 9(3) 


1.669 925 


2.309 81(1) 


2.309 811 


2.5 


2.132 4(8) 


2.132 544 


1.631 83(6) 


1.631 818 


2.452 048(9) 


2.452 048 


3.0 


2.189 29(9) 


2.189 298 


1.622 6(6) 


1.622 556 


2.583 09(7) 


2.583 082 


3.5 


2.244 59(7) 


2.244 588 


1.628 84(6) 


1.628 818 


2.703 9(5) 


2.703 831 


4.0 


2.297 6(3) 


2.297 688 


1.643 9(7) 


1.643 856 


2.815 3(8) 


2.815 465 


4.5 


2.348 3(1) 


2.348 392 


1.664(2) 


1.663 976 


2.919 1(4) 


2.919 098 


5.0 


2.396(1) 


2.396 711 


1.687(7) 


1.687 033 


3.016(3) 


3.015 708 



Since ei vanishes in the ordered phase, the first relation in (3) provides a recursion for the 
components of After normalization one obtains the surface magnetization of the semi- 
infinite system as 

\ -1/2 

CO ] \ 

With an aperiodic system, making use of equations (^ and (H), it may be rewritten as: 

oo 

= [S{\, r)r^'^ , S{\ r)=Y, A-^V'^"^ , no = . (16) 

For the three sequences one obtains a second order surface transition with continuously 
varying exponents Ps{f) (see appendix A for the PF sequence, appendix B for the TF sequence 
and reference for the PD sequence). This nonuniversal behaviour is linked to the value of 
u) = which, together with v~l for the 2d Ising model, leads to a marginal behaviour for the 
three sequences. One may notice that v has to keep its unperturbed value in order to have 
varying exponents. Otherwise the marginality condition (/) = in equation ( p^ ) would no longer 
be satisfied in the perturbed system. Thus, as indicated in the appendices, (3s is equal to the 
scaling dimension of the surface magnetization Xms with: 

^ ln(ri/3+r-i/3) _ln(l+r-i) _ln(2 + r) 

2h^2 ^ 2 1n2 ^ 2 In 3 ^ 

(17) 

on the left surface. The variations with r are shown in figure 2. 

These exponents are the extreme anisotropic limits of those obtained analytically on classical 
2d systems with a finite value of the anisotropy ratio Ki / K2 . Taking the quantum chain 
limit does not change qualitatively the critical behaviour. 
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The magnetic scaling dimension Xms on the right surface is obtained by exchanging perturbed 
and unperturbed couphngs, so that: 

^™.=.„. (PD), = ^ (PF), ^^s = ^^ (TF). (18) 

The two dimensions are the same for the PD sequence which, apart the last digit, is symmetric. 

Let us now consider the surface energy es, with scaling dimension Xes, which is given by 
the ground-state expectation value of erf, i.e. Cg = (0|crf |0). Since Cs contains a non-singular 
part, it is generally difficult to extract its scaling behaviour. One may consider instead the 
connected surface energy-energy correlation function G|(r) = (erf (O)crf (r)) — with scaling 
dimension 2xes- Working in the diagonal basis, one obtains: 

Glir) ^ ^|(^|af|0)|2e-(^-^°) 
j:>0 

a I3>a 

Only pairs of excitations are involved in the second expression. The reason for this is that 
erf only couples the ground state to excited states with two fermions. The first non-vanishing 
matrix element of erf scales like the surface energy. It involves the ground state and the lowest 
two-fcrmions excited state l^} ='7i?72|0) ^'^id, according to (19), may be written as: 

|(e|f7f|0)| = (e2-ei)<?i(l)<?2(l). (20) 

At the critical point, on a chain with length L, this quantity behaves as L~^== . The exponents 
given in table I were deduced from sequence extrapolations using the BST algorithm with 
chain sizes L = 2^, 2'*,- ■ •, 2^6 for the PD and PF sequences and L = 3, 3^,- • •, 3^" for the TF 
sequence. The extrapolated values are consistent with the scaling relation 

For the right surface, x has to be changed into x. 

3.2 Low-lying fermion excitations. — Assuming that the surface component of the sec- 
ond eigenvector in equation (|2^) has the same scaling behaviour as the first one, i.e. that 
$i(l)<I>2(l) ~ one is led to conjecture that low-lying excitations scale like L^^ with, 

according to (|lj) 

Z — Xms ^ras • (^^) 

This has been verified numerically for the six lowest excitations. The results from sequence 
extrapolations, with the same sizes as in the last section, are given in table II for the first 
excitation of the PF sequence. 

Such a behaviour is consistent with a strong anisotropy, i.e. with different correlation length 
exponents, v in the space direction (perpendicular to the layers in the 2d classical system) 
and v\\ = zv in the time direction (along the layers) [Q. The dependence of the anisotropy 
exponent z on r is shown in figure ^ for the three sequences. 

One expects the following behaviour for the excitations under a length rescaling by a factor h: 



(23) 



Fig. 3. — Variation of the anisotropy exponent z with the aperiodic modulation factor r. 



where ts is the surface thermal scaling field, conjugate to the surface energy, with dimension yta- 
Due to the discrete scale invariance of the sequences, the scaling factor has to take values of 
the form m" (m, n integers) where to = 4, 2, 3 for the PD, PF and TF sequences, respectively. 

At the critical point, as shown in table II, the first derivatives with respect to the bulk and 
surface scahng fields scale as: 

These results are consistent with ( p3| ) if as found in the last section, and 

yts — ^ ^es — 2x7ns , (^^) 



in agreement with equation (21). 



Fig. 4. — Test of the anisotropic finite-size-scaling behaviour of the free energy density for the three- 
folding sequence with m — 3 and chain sizes L = 3^ to 3^". 



Table III. — SpeciBc heat exponent for the three marginal sequences studied as function of the 
modulation factor r. For each sequence a is the extrapolated numerical value. The next column gives 
the conjectured analytical value 1 — z. 





period-doubling 


paper 


-folding 


three-folding 


r 


a 


1- z 


a 


1 - z 


a 


1- z 


0.25 


-0.151(5) 


-0.148 844 


-0.29(5) 


-0.321 928 


-0.181(4) 


-0.184 535 


0.5 


-0.039(4) 


-0.038 170 


-0.08(1) 


-0.084 963 


-0.046(3) 


-0.047 952 


0.75 


-0.009(5) 


-0.006 623 


-0.011(5) 


-0.014 874 


-0.009(2) 


-0.008 351 


1.5 


-0.015(3) 


-0.013 137 


-0.028(2) 


-0.029 447 


-0.015(4) 


-0.016 552 


2.0 


-0.045(8) 


-0.038 170 


-0.081(4) 


-0.084 963 


-0.050(4) 


-0.047 952 


2.5 


-0.069(6) 


-0.066 272 


-0.152(7) 


-0.146 391 


-0.086(6) 


-0.082 978 


3.0 


-0.11(3) 


-0.094 649 


-0.207(2) 


-0.207 519 


-0.117(8) 


-0.118 109 


3.5 


-0.13(1) 


-0.122 294 


-0.26(1) 


-0.266 247 


-0.157(7) 


-0.152 102 


4.0 


-0.149(3) 


-0.148 844 


-0.31(1) 


-0.321 928 


-0.186(3) 


-0.184 535 


4.5 


-0.172(4) 


-0.174 196 


-0.36(3) 


-0.374 469 


-0.217(2) 


-0.215 310 


5.0 


-0.193(7) 


-0.198 356 


-0.45(4) 


-0.423 998 


-0.244(3) 


-0.244 464 



3.3 Free energy density and specific heat. — The free energy per row of the system 
is the ground-state energy of the quantum chain which, according to equation (|^), is given by: 

^o = -i$]ea. (26) 

a 



□ 



□ 



Fig. 5. — Scaling functions for the surface spin-spin correlation function. For the period-doubling 
(left) and paper- folding (middle) sequences, r = 5, 4, 3, 2, 1 from top to bottom and L — 2^ (squares), 
2* (triangles) and 2^" (circles). For the three-folding sequence (right), r — 1,2,3,4,5 from top to 
bottom and L = 3* (squares), 3^ (triangles) and (circles). 

At the critical point, with Jl = /(O, 0, 0, L~^) a hnear variation of L {fmL — fh) with 
is expected. In the same way, taking the first derivatives with respect to the bulk and sur- 
face coupHngs, one expects linear variations for L {dfmL/dt— dfL/dt) versus L^~^ and for 
L{dfrnL/dts — dfL/dts) versus l-^-'^'^^^. This has been verified numerically for the three 
sequences. The results for the TF sequence with m = 3 are shown in figure 4. 



Fig. 6. — Scaling functions for the surface energy-energy correlation function. The values of r and L 
are the same as in figure 5. 



The second derivatives of / with respect to t at the critical point has a leading contribution 
from its singular part which behaves as giving the critical exponent of the bulk specific 

heat 



The regular part of the free energy density contributes a constant plus an term which can 
be neglected as long as z < 2. The numerical estimates of a, obtained by extrapolating the 
slopes of In (i9^/mL/9t^— d'^fL/dt^) versus InL, are shown in table III with m = 4 and L^2^ 
to 2^^ for the PD and PF sequences, m — 3 and L==3 to 3^" for the TF sequence. 

3.4 Surface correlation functions. — Once the excitation spectrum and the corre- 
sponding eigenvectors are known, one may calculate the surface spin-spin and energy-energy 
correlation functions as given in equations ( p^ and (19). This has been done at the critical 
point, on finite chains with size L, where the correlation functions scale like: 



a = I — z . 



(29) 



G^(r,L) = b-^'^'^Gl{b-'T,b-^L) , 



(30) 



which gives: 



G^(r,i) 



Gt{r,L) 




= r-2-™=/^ [/.(0) + -^/;(0) + ... 



(31) 



Figures 5 and 6 show the scaling functions for the surface spin-spin and energy-energy 
correlation functions, respectively. Table IV gives the decay exponents deduced from the 
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Table IV. — Decay exponent for the spin-spin (r;"^) and energy-energy (rjj^) surface correlation func- 
tions deduced from sequence extrapolations on the slopes of log-log plots at the critical point. Each 
second column gives the conjectured analytical values. 





period-doubling 


paper 


-folding 


three-folding 


r 




'^^ms / 




2^ms / ^ 




'^■^ms / ^ 


0.5 


1.04(5) 


1 


1.459(2) 


1.460 845 


0.76(5) 


0.795 880 


1 


0.998(7) 


1 


0.98(2) 


1 


0.98(5) 


1 


2.0 


1.02(3) 


1 


0.55(2) 


0.539 155 


1.22(3) 


1.204 120 


3.0 


1.05(7) 


1 


0.36(2) 


0.343 711 


1.30(2) 


1.310 225 


4.0 


1.04(6) 


1 


0.249(7) 


0.243 529 


1.34(4) 


1.376 852 


5.0 


1.05(5) 


1 


0.188(7) 


0.184 715 


1.47(5) 


1.423 298 


r 


vl 






'^•^es / ^ 




2Xe5 / Z 


0.5 


4.01(4) 


4 


4.92(4) 


4.921 691 


3.54(5) 


3.591 760 


1 


3.997(5) 


4 


4.06(7) 


4 


4.01(2) 


4 


2.0 


4.01(2) 


4 


3.13(6) 


3.078 309 


4.42(2) 


4.408 240 


3.0 


4.03(5) 


4 


2.686(5) 


2.687 422 


4.613(9) 


4.620 449 


4.0 


4.02(4) 


4 


2.47(4) 


2.487 058 


4.75(2) 


4.753 704 


5.0 


4.04(6) 


4 


2.34(5) 


2.369 431 


4.88(5) 


4.846 596 



slopes of log-log plots using sequence extrapolations for t = 10 to 100 and the maximum 
chain sizes, L = 2^° for the PD and PF sequences, L = 3^ for the TF sequence. With the 
symmetric PD sequence, Xms = Xms and the decay exponents remain constant in agreement 
with equations (^), (^ and (31). They keep their unperturbed values V^ — ^ and 'yji =4. 

The scaling functions in figures 5 and 6 are clearly oscillating for the PD and TF sequences. 
Such a behaviour is linked with the discrete scale invariance. Renormalization group arguments 
show that critical amplitudes may be log-periodic in such systems (sec for instance [0). When 
the scaling factor is allowed to take values of the form h ~ m" , the amplitude is a periodic 
function with period 1 in the variable In^/lnm where 9 is some scaling field. 

Here, one expects the correlation function to be modulated by an amplitude A which, due 
to the anisotropic scaling, is a periodic function with period 1 in the variable lnT/(z Inm). 
Multiplying both sides of (31) by an oscillating amplitude A, it may be transformed on the 
right into: 



A 



MtL- 



; Ini 



A[^] (32) 
~ In j~ 



when L — m", thus explaining the good data collapse obtained for the scaling functions. 

One may notice that m is generally equal to the leading eigenvalue of the substitution matrix 
Ai [Q. Although this is verified for the PF and TF sequences, m = Af for the PD sequence. It 
may be verified on (|^) that the behaviour of this sequence is different for odd and even values 
of n. 



3.5 Bulk magnetization. — On a chain with free boundary conditions, the local magneti- 
zation on site I is a matrix element of af which can be evaluated in the diagonal basis, as in 



Fig. 7. — Bulk magnetic exponent deduced from finite-size scaling and sequence extrapolations as a 
function of the modulation amplitude for different values of the ratio 1/ L. The circles for the period- 
doubling sequence are obtained when the unperturbed coupling at l = L/2 is replaced by a perturbed 
one. 



section 3.1, giving: 



mi = ((7|af|0) ^ (0|??i^iSiA2S2' 

Q. 



■ Ai^iBi^iAi\0) , 



(33) 



The local magnetization mi can be calculated, like correlation functions in reference pd[ , by 
means of Wick's theorem. The ground-state expectation value is given by sums of products of 



pair expectation values. Since {0\AjAk\0) -- 
is obtained as the following determinant: 



\BjBk\0) =0 when j ^k, the local magnetization 



mi 



Hi Gil Gi2 
H2 G2 1 G2 2 

Hi Gil G12 



where: 



Gi i-i 
G21-1 



G 



I i-i 



(34) 



G 



(0|r;iA,|0)=<?i(i), 
(0|BfcA,|0) = -^<l'„(j)tf'„(fc) 



(35) 



Fig. 8. — Scaled energy density profile for the three sequences with a modulation factor r — 2 and 5 
and chain sizes L = 2^ (squares), 2* (triangles), 2^" (circles) for the period-doubling and paper-folding 
sequences, L = 3* (squares), 3^ (triangles), 3^ (circles) for the three-folding sequence. 



At the critical point, when L = 77i", the magnetization is expected to behave as: 

mi{L-^) = 6-^" mi/t{bL-^) (36) 

so that keeping the ratio l/L constant, Xm governs the size-dependance of the magnetization, 
mi ~ L"^". 

The bulk magnetization exponent Xm has been obtained through finite-size scaling and se- 
quence extrapolations on chains with sizes 2^, 2*,- ■ •, 2^" for the PD and PF sequences and 
3,3^,- • -,3^ for the TF sequence. The results are shown in figure 7. The exponent is continu- 
ously varying with r due to the marginal aperiodicity. It depends also on the ratio l/L, which 
means that the scaling relation (|36| ) is only valid locally with the local value of Xm- 

With the PD sequence, for example, different exponents are obtained at l — L/2 on one hand 
and l — L/A or 3L/A on the other hand. This can be understood by looking at the sequence 
near these sites for long enough chains with sizes of the form L^". In the first case, one finds 
the following sequence, 1011101010111010, whereas the environment is the same but 
is replaced by 1 in the other case. This is enough to change the local critical behaviour because 
a thermal line defect is a marginal perturbation when i/ = 1 and d = 2 . If one changes 
the (underlined) unperturbed coupling A at L/2 by a perturbed one Ar, the local exponent 
becomes the same as for L/4 and 3i/4. The same type of behaviour is obtained with the two 
other sequences (see figure 7). This inhomogeneous critical behaviour suggests a multifractal 
statistics as for electronic states on the same structures |53) . 

3.6 Bulk energy density and its profile. — We now consider finite chains with length L 
and periodic boundary conditions, i.e. with <l>a{L + l) = — Q(l>a{l) and Xo^a{0) =^ —Q^L'I^a{L) 
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Table V. — Bulk energy exponent for the three marginal sequences studied as function of the mod- 
ulation factor r. For each sequence Xe is the numerical value deduced from sequence extrapolations 
whereas the next column gives the conjectured analytical value. The figure in brackets gives the 
uncertainty on the last digit. 





period-doubling 


paper-folding 


three-folding 


r 




z 




z 




z 


0.5 


1.03(4) 


1.038 170 


1.084 9(7) 


1.084 963 


1.048(5) 


1.047 952 


1. 


0.999 99(2) 


1 


0.999 98(3) 


1 


1.000(1) 


1 


1.5 


1.02(2) 


1.013 137 


1.029 4(2) 


1.029 447 


1.016 5(6) 


1.016 552 


2.0 


1.06(3) 


1.038 170 


1.084 95(3) 


1.084 963 


1.047 8(6) 


1.047 952 


2.5 


1.09(4) 


1.066 272 


1.146 38(5) 


1.146 391 


1.083(1) 


1.082 978 


3.0 


1.12(5) 


1.094 649 


1.207 51(7) 


1.207 519 


1.118(2) 


1.118 109 


3.5 


1.15(5) 


1.122 294 


1.266 24(2) 


1.266 247 


1.152(1) 


1.152 102 


4.0 


1.18(6) 


1.148 844 


1.321 92(2) 


1.321 928 


1.184(1) 


1.184 535 


4.5 


1.20(6) 


1.174 196 


1.374 5(6) 


1.374 469 


1.215(1) 


1.215 310 


5.0 


1.22(6) 


1.198 356 


1.424 0(9) 


1.423 998 


1.244(2) 


1.244 464 



where the conserved charge Q = Ofe=i '^k commutes with the Hamiltonian. Thus one may 
proceed by replacing the charge operator through one of its eigenvalues, Q = ±l provided the 
levels with the appropriate charge are retained in each sector (see for instance reference fs^). 

It has been verified in both sectors that low-lying single-particle excitations scale as , like 
in the system with free boundaries. Furthermore, the free energy density can be analysed as 
in section 3.3 and is consistent with anisotropic scaling. One may notice that there remains 
a somewhat unexpected contribution to the free energy density which does not appear 
in homogeneous systems. It seems that closing the aperiodic sequence onto itself is similar to 
introducing a defect which contributes a size-independent term to the free energy. 

We shall discuss more extensively the scaling of the bulk energy density which, on site 
follows from the matrix element 

ei = (e|af|0) = \Mm2{l)-Ml)^i{l)\ , (37) 

where the ground-state |0) and the two-fermion eigenstate \e) = ?y|?7||0) both belong to the 
even sector. The following behaviour is expected: 

ei{t,L-^)^b~'ei/t{b^^''t,bL~^) (38) 

with b = m'^ and i/ = l as already mentioned. 

According to ( ^Sf ) , at the critical point the scaled energy density profile ei is a function 
of l/L only. This is verified in figure 8 for the three sequences. A good data collapse is 
obtained as long as r is not too large and there is no sign of log-periodic modulation in this 
case. One may notice that the scaling assumption ( p8| ) supposes that the energy density 
keeps the same scaling dimension z everywhere on the periodic chain, i.e. it does not display 
the inhomogeneous behaviour observed for the magnetization. Only the amplitude of is 
varying with l/L. The bulk energy exponent deduced from finite-size scaling on the first site 
is compared to the conjectured analytical result for the three sequences in table V. 
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4. Conclusion. 

We have presented the resuhs of a systematic study of the influence of marginal aperiodic 
perturbations on the critical behaviour of the Ising quantum chain. For the three marginal 
sequences considered, the surface magnetization has been obtained exactly. The surface mag- 
netic exponents are continuously varying with the modulation amplitude and different values 
Xms and Xms are obtained on the two surfaces when the sequence is asymmetric. 

The critical behaviour of the excitation spectrum, the free energy density, the specific heat 
and the surface correlation functions on systems with free boundary conditions have been 
studied numerically through finite-size scaling. The results are in agreement with anisotropic 
scaling with an anisotropy exponent z which is simply given by the sum of the two magnetic 
surface exponents. The thermal surface and bulk exponents have been also related to Xms and 
Xms- As a consequence of the discrete scale invariance, the scaling functions of the surface 
correlation functions display log-periodic modulations for the PD and TF sequences. 

The bulk energy proflle and the bulk magnetization obtained on systems with periodic bound- 
ary conditions also display anisotropic scaling. The bulk magnetic exponent is continuously 
varying and may take different values around sites corresponding to simple fractions of the 
chain length. Although the bulk energy exponent is continuously varying too, this variation is 
a simple consequence of anisotropic scaling: Xe — z, which is varying. The correlation length ex- 
ponent itself keeps its unperturbed value. This is consistent with the occurence of continuously 
varying exponents, since the value i^=l is necessary to satisfy the marginality condition. 

One may notice that the anisotropy exponent z is always > 1 for the three sequences studied. 
As a consequence the specific heat exponent a = 1 — z < and the Ising logarithmic singularity 
is washed out. 

The anisotropic critical behaviour associated with these marginal systems breaks the con- 
formal invariance of the critical Ising model. It follows that the gap-exponent relations are no 
longer valid. 

Much remains to be done about these systems. One needs some proofs for the different 
conjectures about the values of the critical exponents. In particular, an analytical calculation of 
the anisotropy exponent would be useful. In this paper, we did not examine systematically the 
log-periodic oscillations of the critical amplitudes. The possibility for a multifractal behaviour 
of the bulk magnetization has also to be explored. Finally, it would be interesting to consider 
in detail the case of relevant aperiodic perturbations. 
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Appendix A. 
Paper-folding sequence. 

The paper-folding sequence jii) results from the recurrent folding of a sheet of paper, right 
over left, as shown for the first few steps in figure 9. After unfolding, one obtains a succession 
of up- or down-folds. Associating 1 with an up-fold and with a down-fold leads, after four 
foldings, to the following sequences: 

110110011100100 (A.l) 



Fig. 9. — Folding n times the right part of a sheet of paper onto the left gives, after unfolding, a 
succession of up- or down-folds which follow the paper-folding sequence. 



The sequence on the right of the central fold is the mirror image of the sequence on the left, 
with each digit fk replaced by its complement 1 — fk- As a consequence, the asymptotic density 
is Poa — ^ and equation (^) gives the critical coupling Xc = r^^^^. 

The same sequence can be generated using the two-digit substitutions <S(00) = 1000, 5(01) = 
1001, 5(10) = 1100, 5(11) = 1101. Starting with 11, the resulting sequence differs from the 
previous one only through its last supplementary digit, 1, which does not affect the critical 
behaviour. 

According to (p]), the leading eigenvalues of the substitution matrix, Ai =2, A2 = 1, lead to 
a vanishing wandering exponent, w = 0, i.e. the corresponding layered perturbation in the 2d 
Ising model is marginal. 

Underlined even terms in ( AA ) reproduce the sequence whereas odd terms are alternatively 
1 or so that we have: 

f2k = fk , /2fc+i ^ 2 I-"*" (^1)''] ' 

p 3 + (-i)p p 1 - i-iy 

The sum in ( p^ can be rewritten as 

OG 00 

p=0 p=0 

Separating odd and even values of p in the two sums and using the recursion relations in (A. 2), 
we obtain the following set of coupled functional equations: 

5odd(A,r) = A-2r-iS'„dd(AVi/2,r) + A-V-25evcn(A2ri/2,r) 

5even(A,r) =r-i5odd(AVi/2,r) + 5cvcn(AVi/2,r) . (A.4) 

These relations can be written in matrix form giving Vq = rip=o 2!^ YLk where: 



2" 



638 



JOURNAL DE PHYSIQUE I 



N°5 



In the ordered phase A > Ac and hnip^oo }Lz 



5(A,r) = (l 1) 



.P=0 



so that, according to (A. 4), 



(A.6) 



which, together with equation (|16D, formally gives the surface magnetization. 

The critical behaviour can be obtained using the scaling method of reference [^. Let S{z) 
be the series expansion of S{X,r) in powers of z= (Ac/A)^. According to equation (|l6|), near 
the critical point S{z) behaves as (1 — z)~'^^' . Then it may be shown that the truncated series 
Sl{z), containing the first L terms in S{z), diverges as L"^^' at the critical point z = l when 
L ^ oo. According to finite-size scaling, at the critical point of a system with size L, iris 
vanishes as L"^""" so that S'l(I) also behaves as L^^™=. As a consequence (3s = Xms and ^ = 1 
like in the pure system. 

The form of the matrix in ( |A.5[ ) is such that, keeping the first n factors in the infinite product 
of equation (A.6), gives 



SUz) = (1 1) 



n 

.p=0 



2" 



Z 

A? 



Xlz^^ 
1 



(A.7) 



with L = 2". Since at z = 1 the matrix becomes p-independent, the form of the leading 
contribution to the truncated series is easily obtained as 

5l=2"(1)-^(2")2^=^(1 + A2)", (A.8) 

where the last term is the nth power of the leading eigenvalue of T . Finally, the surface 
magnetization exponent reads 

ln(l + r-i) 



2 In 2 



(A.9) 



Appendix B. 
Three-folding sequence. 

The three-folding sequence Q is generated through the substitutions 5(0) = 010, 5(1) = 011. 
Starting with 0, after four steps one obtains: 

010011010010011011010011010 (B.l) 

The sequence repeats itself if one keeps every third (underlined) term, which gives: 

fsk = fk , fsk+l — , /3fe+2 = 1 • 

«3fe ^rik + k , ri'ik+i = nfe + fc , ns,k+2 = Uk + k + I . (B.2) 

Iterating 713^ gives n^p — ^{V — 1) so that, asymptotically, poo = \ and Ac = The 
eigenvalues of the substitution matrix, Ai = 3, A2 = 1, lead to the wandering exponent w = 0, 
i.e. marginal behaviour in the 2d Ising model. 

Like in (A. 3) the sum ^(A, r) may be splitted into three sums which, making use of (B.2), 



lead to the recursion: 



S{\r) 



(l + A-2+r-2A-4)5(AV, r) 



n(i+A 

p=0 



-2 
p 



(B.3) 
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where Xp = Ac {X/Xc)^" . Using the same scaling method and notations as in appendix A, the 
first n factors in the infinite product contain the first L = 3" terms in the truncated series, 



p=0 

so that, from its value at the critical point z = l, one deduces 

ln(2+r) 



2 ln3 



(B.4) 



(B.5) 
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